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Abstract 

We study the competition of quark-antiquark and diquark condensates under 
the influence of an external chromomagnetic field modelling the gluon conden- 
sate and in dependence on the chemical potential and temperature. As our 
results indicate, an external chromomagnetic field might produce remarkable 
qualitative changes in the picture of the color superconducting (CSC) phase 
formation. This concerns, in particular, the possibility of a transition to the 
CSC phase and diquark condensation at finite temperature. 
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I. INTRODUCTION 



Low energy (large distance) effects in QCD can only be studied by approxi- 
mate (nonperturbative) methods in the framework of various effective models 
or in terms of lattice calculations. At present time, one of the most popu- 
lar QCD-like effective theories is the well-known Nambu-Jona-Lasinio (NJL) 
model which is a relativistic quantum field theory with four-fermion in- 
teractions. The physics of light mesons (see e.g. [|2] and references therein), 
diquarks j3|,[|] and meson-baryon interactions j5j]- |7j based on dynamical 
chiral symmetry breaking can be effectively described by NJL chiral quark 
models. Moreover, NJL models are widely used in nuclear physics and astro- 
physics (neutron stars) for the investigation of quark matter Q8j , to construct 
alternative models of electroweak interactions Q and in cosmological applica- 
tions pDf] . Moreover, its (2+l)-dimensional analogue serves as a satisfactory 



microscopic theory for several effects in the physics of high-temperature su- 
perconductors 



The NJL model displays the same symmetries as QCD. So it can be suc- 
cessfully used for simulating some of the QCD vacuum properties under the 
influence of external conditions such as temperature T and chemical potential 
ji . The role of such considerations significantly increases especially in the 



cases, where numerical lattice calculations are not admissible in QCD, i.e. at 
nonzero density and in the presence of external electromagnetic fields [jT3| , [H| . 



Recently, it was shown in the framework of a (2+l)-dimensional NJL model 
that an arbitrary small external magnetic field induces the spontaneous chi- 
ral symmetry breaking (xSB) even under conditions, when the interaction 
between fermions is arbitrary weak fpT5fl - Later it was shown that this phe- 
nomenon (called magnetic catalysis effect) has a rather universal character 
and gets its explanation on the basis of the dimensional reduction mechanism 
jn|. (The recent reviews flT7j consider the modern status of the magnetic 
catalysis effect and its applications in different branches of physics.) 

As an effective theory for low energy QCD, the NJL model does not con- 
tain any dynamical gluon fields. Such nonperturbative feature of the real 
QCD vacuum, as the nonzero gluon condensate < F^F afll/ >=< FF > can, 
however, be mimicked in the framework of NJL models with the help of exter- 
nal chromomagnetic fields. In particular, for a QCD-motivated NJL model 
with gluon condensate (i.e. in the presence of an external chromomagnetic 
field) and finite temperature, it was shown that a weak gluon condensate 
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plays a stabilizing role for the behavior of the constituent quark mass, the 
quark condensate, meson masses and coupling constants for varying temper- 
ature Then, in a series of papers, devoted to the NJL model with gluon 
condensate, it was shown that an external chromomagnetic field, similar to 
the ordinary magnetic field, serves as a catalyzing factor in the fermion mass 
generation and dynamical breaking of chiral symmetry as well |H| . The basis 
for this phenomenon is the effective reduction of the space dimensionality in 



the presence of external chromomagnetic fields pO 



There exists the exciting idea proposed more than twenty years ago [^T 



2Bj that at high baryon densities a colored diquark condensate < qq > might 
appear. In analogy with ordinary superconductivity, this effect was called 
color superconductivity (CSC). The CSC phenomenon was investigated in 
the framework of the one-gluon exchange approximation in QCD [Z3j] , where 
the colored Cooper pair formation is predicted selfconsistently at extremely 
high values of the chemical potential fi ^ 10 8 MeV RZ5fl . Unfortunately, such 
baryon densities are not observable in nature and not accessible in experi- 
ments (the typical densities inside the neutron stars or in the future heavy 
ion experiments correspond to /i ~ 500 MeV). The possibility for the exis- 
tency of the CSC phase in the region of moderate densities was proved quite 
recently (see e.g. the papers [ffijfl - [|29|1 as well as the review articles [j3t| and 



references therein). In these papers it was shown on the basis of different 
effective theories for low energy QCD (instanton model, NJL model, etc.) 
that the diquark condensate < qq > can appear already at a rather moder- 
ate baryon density (/i ~ 400 MeV), which can possibly be detected in the 
future experiments on ion-ion collisions. Since quark Cooper pairing occurs 
in the color anti-triplet channel, the nonzero value of < qq > means that, 
apart from the electromagnetic U(l) symmetry, the color SU C (3) should be 
spontaneously broken down inside the CSC phase as well. In the framework 
of NJL models the CSC phase formation has generally be considered as a 
dynamical competition between diquark < qq > and usual quark- antiquark 
condensation < qq >. However, the real QCD vacuum is characterized in 
addition by the appearence of a gluon condensate < FF > as well, which 
might change the generally accepted conditions for the CSC observation. In 
particular, one would expect that, similarly to the case of quark- antiquark 
condensation, the process of diquark condensation might be induced by ex- 
ternal chromomagnetic fields. For a (2+l)-dimensional quark model, this was 
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recently demonstrated in [[Ttj. There, a SU (2) L x SU (2) R chirally symmetric 
(2+l)-dimensional NJL model with three colored quarks of two flavors was 
considered at zero T, fi. It was shown that in this case for arbitrary fixed 
values of coupling constants there exists a critical value of the external chro- 
momagnetic field at which a CSC second order phase transition is induced in 
the system. Q Since the two-flavored QCD3 and the considered NJL model 
are not in the same universality class of theories (QCD 3 with Nj = 2 has 
a higher flavor symmetry S77(4)), the obtained results are intrinsic to real 
QCD4 rather than to QCD3. Indeed, our recent investigations on the basis of 
a (3+l)-dimensional NJL model [£52|,E53j and \i = show that some types of 
sufficiently strong external chromomagnetic fields may catalyze the diquark 
condensation. 

As argued above, CSC might occur inside neutron stars and possibly be- 
come observable in ion-ion collisions, i.e. at nonzero baryon densities. Taking 
into account the fact that at finite chemical potential the magnetic gener- 
ation of dynamical x>3B qualitatively differs from the /i = case |T3j| , one 
might expect analogous effects for CSC, too. By this reason, the investiga- 
tion of the chromomagnetic generation of CSC under the influence of a finite 
chemical potential (finite particle density) is a very interesting and actual 
physical problem. 

The aim of the present paper is to study the influence of external con- 
ditions such as chemical potential, temperature and especially of the gluon 
condensate (as modelled by external color gauge fields) on the phase struc- 
ture of quark matter with particular emphasize of its CSC phase. To this 
end, we shall extend our earlier analysis of the chromomagnetic generation 
of CSC at fi = |3T]]- |33[ to the case of an (3+l)-dimensional NJL type 



model with finite chromomagnetic field, temperature and chemical potential 
presenting a generalization of the free field model of 



The paper is organized as follows. In Sections II and III the extended NJL 
model under consideration is presented, and its effective potential (= ther- 
modynamic potential) at nonzero external chromomagnetic field, chemical 
potential and temperature is obtained in the one-loop approximation. This 
quantity contains all the necessary informations about the quark and diquark 



Strictly speaking, the CSC is induced by those components of external chromomagnetic fields, 
which can stay massless inside the CSC phase. 
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condensates of the theory. In the following Section IV the phase structure of 
the model is discussed on the basis of numerical investigations of the global 
minimum point of the effective potential. As our main result, it is shown 
that the external chromomagnetic field can induce the transition to the CSC 
phase and diquark condensation even at finite temperature. Thereby, the 
characteristics of the CSC phase can significantly change in dependence on 
the strength of the chromomagnetic field. Finally, Section V contains a sum- 
mary and discussion of the results. Some details of the effective potential 
calculation are relegated to an Appendix. 

II. THE MODEL 

Let us first give several (very approximative) arguments motivating the 
chosen structure of our QCD-motivated extended NJL model introduced be- 
low. For this aim, consider two-flavor QCD with nonzero chemical potential 
and color group SU C (N C ) and decompose the gluon field A®(x) into a con- 
densate background ("external") field A% (x) and the quantum fluctuation 
a"(x) around it, i.e. A%(x) =A®(x)-\-a®(x). By integrating in the generating 
functional of QCD over the quantum field a%(x) and further "approximat- 
ing" the nonperturbative gluon propagator by a function, one arrives at 
an effective local chiral four-quark interaction of the NJL type describing 
low energy hadron physics in the presence of a gluon condensate. Finally, by 
performing a Fierz transformation of the interaction term, one obtains a four- 
fermionic model with (^)-and (gg)-interactions and an external condensate 
field A a (x) of the color group SU C (N C ) given by the following Lagrangian Q 

L = q[Y(id v + gAl{x)^) + py°\ q + ^[(qqf + (qi^rq) 2 ] + 

+ ^[igce(iAL)7 5 9][i?e(iAL)7 5 9c]. (1) 

It is necessary to note that in order to obtain realistic estimates for masses 
of vector/axial- vector mesons and diquarks in extended NJL-type of models 



2 The most general four-fermion interaction would include additional vector and axial-vector (qq) 
as well as pseudo-scalar, vector and axial-vector-like (qq) -interactions. For our goal of studying 
the effect of chromomagnetic catalysis for the competition of quark and diquark condensates, the 
interaction structure of (1) is, however, sufficiently general. 
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[J3j], we have to allow for independent coupling constants G\, G^-, rather than 
to consider them related by a Fierz transformation of a current-current in- 
teraction via gluon exchange. Clearly, such a procedure does not spoil chiral 
symmetry. 

In (UD g denotes the gluon coupling constant, \x is the quark chemical 
potential, q c = C(f, q c = q t C are charge-conjugated spinors, and C = i^ 2 ^ 
is the charge conjugation matrix (t denotes the transposition operation). In 
what follows we assume N c = 3 and replace the antisymmetric color matrices 
A^ s (with a factor i) by the antisymmetric e b operator. Moreover, summation 
over repeated color indices a = l,...,8;6=l,2,3 and Lorentz indices v = 
0, 1, 2, 3 is implied. The quark field q = qi a is a flavor doublet and color triplet 
as well as a four-component Dirac spinor, where % = 1, 2; a = 1, 2, 3. (Latin 
and Greek indices refer to flavor and color indices, respectively; spinor indices 
are omitted.) Furthermore, we use the notations A a /2 for the generators 
of the color SU C (3) group appearing in the covariant derivative as well as 
r = (r 1 , t 2 , r 3 ) for Pauli matrices in the flavor space; {e) lk — e tk , (e b ) af3 — e af3b 
are totally antisymmetric tensors in the flavor and color spaces, respectively. 
Clearly, the Lagrangian (jlj) is invariant under the chiral SU(2)l X SU(2)r 
and color SU C (3) groups. 

Next, let us for a moment suppose that in ([I]) A a (x) is an arbitrary clas- 
sical gauge field of the color group SU C (3). (The following investigations do 
not require the explicit inclusion of the gauge field part of the Lagrangian). 
The detailed structure of A a (x) corresponding to a constant chromomagnetic 
gluon condensate will be given below. 

The linearized version of the model (|TJ) with auxiliary bosonic fields has 
the following form 

L = q[Y(id„ + gAt(x)^) + ^°]q - q(a + i^)q - ^-(a 2 + t? 2 ) - 

- J- A* b A b - ^[iq'Cee^q] - A^iqee^Cq 4 } . (2) 
G2 

The Lagrangians (|I]) and (0) are equivalent, as can be seen by using the 
equations of motion for bosonic fields, from which it follows that 

A b ~ iq t Cee b j 5 q, o ~ qq, tt ~ iqj 5 fq. (3) 

Clearly, the a and 7? fields are color singlets. Besides, the (bosonic) di- 
quark field A b is a color antitriplet and a (isoscalar) singlet under the chiral 
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SU(2)l x SU(2)r group. Note further that the a, A 6 , are scalars, but the tt 
are pseudo-scalar fields. Hence, if a ^ 0, then chiral symmetry of the model 
is spontaneously broken, whereas A b ^ indicates the dynamical breaking 
of both the color and electromagnetic symmetries of the theory. 

In the one-loop approximation, the effective action for the boson fields 
which is invariant under the chiral (flavor) as well as color and Lorentz groups 
is expressed through the path integral over quark fields: 

exp(zSeff(<7,7r, A 6 , A*\A a )) = N' f[dq}[dq] exp{i f Ld A x), 



where 



S eS (a^,A t} ,A*\A a ) = -N c l d*x 



g 2 + k 2 
2G X 



b A *b 



+ 



A b A 



+ S, 



(4) 



N' is a normalization constant. The quark contribution to the partition 
function is here given by: 

Z q = exp(iS) = N' f [dq] [dq] exp(i f [qVq + qMq 4 + q t Mq)d A x). (5) 



In (5) we have used the following notations 



V = D + l»gAl(x)'—- D = iYd» -a- i^r + py°, 
M = -iA* b C£€ b -? 5 , M = -iA b £e b -? 5 C, 



(6) 



where D is the Dirac operator in the coordinate, spinor and flavor spaces, 
whereas V, M. and M. are in addition operators in the color space, too. Let us 
next assume that in the ground state of our model (A 1 ) = (A 2 ) = (7?) = and 
(a), (A 3 ) 7^ 0. Q Obviously, the residual symmetry group of such a vacuum 
is SU C (2) whose generators are the first three generators of the initial SU C (3). 
Now suppose that in this frame the constant external chromomagnetic field, 
simulating the presence of a gluon condensate (FF) = 2H 2 , has the following 
form H a = (H 1 , H 2 , H 3 , 0, . . . , 0). Furthermore, due to the residual SU C (2) 
invariance of the vacuum, one can put H 1 = H 2 = and H 3 — H. 

Some remarks about the structure of the external chromomagnetic fields 
A®(x) used in ([Ij) are needed. From this moment on, we assume A®(x) in 



3 If (if) 7^ then one would have spontaneous breaking of parity. For strong interactions parity is, 
however, a conserved quantum number, justifying the assumption (if) = 0. 



6 



such a form that the only nonvanishing components of the corresponding field 
strength tensor F" are Ff 2 = — F 2 3 : = H = const. The above homogeneous 
chromomagnetic field can be generated by the following vector-potential 

Al(x) = {0,0,Hx\0); A a v {x) = (a ^ 3), (7) 

which defines the well known Matinyan-Savvidy model of the gluon conden- 
sate in QCD (g . 

In QCD the physical vacuum may be interpreted as a region splitted 
into an infinite number of domains with macroscopic extension j35j]. Inside 



each such domain there can be excited a homogeneous background chro- 
momagnetic field, which generates a nonzero gluon condensate (FF) ^ 0. 
(Averaging over all domains results in a zero background chromomagnetic 
field, hence color as well as Lorentz symmetries are not broken.) ^\ 

In order to find nonvanishing condensates (a) and (A 3 ), we should cal- 
culate the effective potential, whose global minimum point provides us with 
these quantities. Suppose that (apart from the external vector-potential 
A®(x) (0)) all boson fields in S e & (H) do not depend on space-time. In this 
case, by definition, S c s = —V e s J d 4 x, where 



Due to our assumption on the vacuum structure, we put A 1 ' 2 = 0, as well 
as 7T = 0. Then, taking into account the form of the vector-potential (7), one 
can easily see that the functional integral for S in (5) is factorized 

Z q = exp(iS(<7, A)) = N' J [dq 3 ] [dq 3 ] exp(i J q 3 Dq 3 d A x) x (9) 
x J [dQ] [dQ] exp(* / [QVQ + QMQ 1 + Q t MQ]d i x), (10) 

where A = A 3 , ^3 is the quark field of color 3 and Q = (q\, qz) 1 is the doublet, 
composed from quark fields of the colors 1,2. Moreover, D = D\j?=Q (D is 
presented in (6)) and 

T> = D + ^gAl(x)^; M = -iA*Cee^ 5 , M = -lAee^C. (11) 
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Strictly speaking, our following calculations refer to some given macroscopic domain. The ob- 
tained results turn out to depend on color and rotational (Lorentz) invariant quantities only, and 
are independent on the concrete domain. 



7 



In ( |TI1) (J3, 1 are matrices in the two-dimensional color subspace, correspond- 
ing to the SU C {2) group: 



^3 



10 „_ / 1 

0-1 ' *~ -10 



Clearly, the integration over q% in (g) yields det D. 

Defining ty* = (Q*, Q) and introducing the matrix-valued operator 



Z = 



2M , -V 1 
V , 2M 



the gaussian integral over Q and Q in ([T0|) can be rewritten in compact 
matrix notation and be evaluated as 

[cM]eif* tzm4x = detV 2 Z. (12) 
Then, by using in (|T2f) the general formula 

det (I ' a) = det [~^ 5 + BAB- 1 A] = det[AA - AbA 1 B], 
one obtains the result: 

exp(iS(a, A)) = iVMet(Z))det 1//2 [4M M + MV l M~ x V\ = 
= Ar'det[(ia-(7 + /i7 )] • 

•det 1/2 [4|A| 2 + (-id - u + M7° - <?i 3 — )(*d - a + /i 7 ° + #i 3 — )] . (13) 

Recall that the operator under the first det-symbol in ( |I3|) acts only in the 
flavor, coordinate and spinor spaces, whereas the operator under the second 
det-symbol acts in the two-dimensional color subspace, too. 

III. THE GENERAL CASE ///0,T^0,iI^0 
A. The effective potential 

First of all, let us calculate the effective action from ([T3|) at zero temper- 
ature T. It is convenient to rewrite the second determinant in (|T31 ) in the 
form 



S 



det[4MM + MV t (M)~ 1 T>] = 
= det [4 1 A| 2 + /i 2 - pi + a 2 - (7V) 2 - 2/i 1 °(a + 7 V)] (14) 

where the /-momentum space representation and 7V = 7fc(idfc + gA\ao,/2), 
(k = 1,2,3) have been used. Similarly to quantum electrodynamics, 
it is easily seen that the operator H = 7°(cr + 7V) is the Hamilto- 
nian for quarks with color indices a = 1,2 and flavor i = 1,2 in the 
background vector-potential (j7|). Its eigenvalues are where £s n \ = 



° 2 + Pi + gH{n + 1/2) — gH(/2, and corresponding eigenstates are de- 
noted by &f n \p 2ia - The set of quark quantum numbers in the background field 
are defined as follows: {n} = {n = 0, 1,2, ...; — 00 < p% < +00; ( = ±1}, 
i,a = 1,2, —00 < P2 < 00. Each of the eigenvalues for Ti is evidently 

four-fold degenerate with respect to flavor and color quantum numbers i, a. 
It is also degenerate with respect to the quantum number P2, which quasi- 
classically characterizes the charged particle center of orbit position in an 
external uniform magnetic field. Since 

one can easily conclude that in the basis $f n y p2ia the operator in the deter- 
minant (|I3j) is diagonal. Moreover, its diagonal matrix elements are equal 
to 4|A| 2 — pi + (/i ± £{ n }) 2 - Upon multiplying these quantities, one can 
find the determinant from QT3Q. In a similar way, it is possible to calcu- 
late the first determinant from flT3Q. Hence, taking into account the relation 
DetO = exp(TrlnO), and following the standard procedure (see, e.g., HffiQ, 
the following expression for V is obtained from ( omitting an infinite 

a- and A-independent constant): 

V = - S - = iN f j g{ E ± In ((E p ± - pi) + 

+ ln(4|A| 2 - 1? 2 + (£ M ± / i) 2 )}, (15) 

{n},± } 

where {p}o denotes the set of quark quantum numbers for vanishing back- 
ground field ({p}o = { — 00 < pi, p2, Ps < +00 }), and E p = \/p 2 + a 2 . The 
factor Nf in front of the integral in ( |T5|) is the result of summation over flavor 
indices % = 1, . . . , Nf, whereas the degeneracy factor A = gH/ (87r 2 ) is due to 
the integration over the momentum pi and summation over the color indices 
a = 1,2. Moreover, £ Wo = Jd 3 p/{27r) 3 , E {n } fdp 3 / {2tt) E n ,c- 
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In the case of finite temperature T = 1/(3 > the corresponding expres- 
sion for Vt can be obtained from flI5Q by means of the following replacements: 

/^(••O-tfTG-O; p ^iu l = 2 7 nT(l + 1/2); I = 0, ±1, ±2, . . . , 
J Ztt l 

where tui is the Matsubara frequency. Hence, 

l=ao , 

V T = -N f T £ { E In ((E p ± + cuf) + 

Z=-00 {p} ,± 

+ 4 £ ln(4|A| 2 + ^ 2 + (£ M ±/i) 2 )}. (16) 

{n},± 

In order to transform ( |T§) , let us first perform the summation over the 
Matsubara frequencies. It is evident that 

£ ln(u,f + n 2 ) = E / da 2 -^ + £ In L 2 + 1) , (17) 



\2 



where Q stands, according to ([16]), for ± fi) 2 or y4|A| 2 + (s{ n } ± /i) 

i.e. Q > 0. Note that we can neglect the contribution from the last term in 
( |I7|) , since it does not depend on o and A. The first term in ( |T7|) can be 
presented in the following form (see Appendix): 

£ / da 2 — 2 = 21nch(Q/?/2) + const 

l v > w z + a 

= Qf3 + 2 ln(l + e~ 0/3 ) + const . (18) 

Performing the summation over Matsubara frequencies in the second term 
in (|I51), and taking into account the degeneracy of the quark spectrum e n in 



the chromomagnetic field with respect to combination of quantum numbers 
n and £, we can use the following expression for the energy spectrum: e n = 



gHn + p§ + a 2 , where n = 0, 1, 2, ... is the Landau quantum number, and 
— oo < £>3 < oo. Then, summing over the spin quantum number £ = ±1, 
we have to account for the fact that for the ground state with n = only 
one spin projection ( = — 1 is possible. Hence, a factor a n = 2 — 5 n o should 
be included in the final expression. As for the summation over Matsubara 
frequencies in the first term in (|Tfi|) , it is necessary to take into account the 
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fact that the function ((T8j) is even with respect to the variable Q. Finally, we 
thus arrive at the following result for the thermodynamic potential 

V H , T (a,A) = N c 

- 2N f J -0y S (N c - 2)[E P + Tln[(l + e -«^)(l + e"^^)]} - 

- N f A £ ^3« n {#n-rt 2 + 4|A| 2 + V / (£ n + A i) 2 + 4|A| 2 + 
+ 2Tln[(l + e -N^-tf+m* ){ i + e-V( £ "+M) 2 +4|A| 2)] \ 




For convenience, expressions are again written in terms of Nf and AT C even 
though in the following we will be concerned only with Nf = 2 and N c = 3. 

B. Regularization 

First of all, let us subtract from ( |I9l) an infinite constant in order that the 
effective potential obeys the constraint Vflyr(0, 0) = 0. After this subtraction 
the effective potential still remains UV divergent. This divergency could 
evidently be removed by introducing a simple momentum cutoff \p\ < A. 
Instead of doing this, we find it convenient to use another regularization 
procedure. To this end, let us recall that all UV divergent contributions to 
the subtracted potential Vh^t(o', A) — Vflyr(0, 0) are proportional to powers 
of meson and/or diquark fields cr, A. So, one can insert some momentum- 
dependent form factors in front of composite cr-and A-fields in order to 
regularize the UV behaviour of integrals and sums, f] 

It is clear by now that we are going to study the effects of an external 
chromomagnetic condensate field in the framework of the NJL-type model 



5 A suitable physical motivation for such form factors follows in the framework of nonlocal NJL type 
models based on the one-gluon-exchange approximation to QCD with nontrivial gluon propagator. 
In particular, in ref. [37] it was shown that the arising exponential form factors for composite 
mesons, as obtained from the solution of the (nonlocal) Bete-Salpeter(BS)-equation, make the 
quark loop expansion including meson (diquark) vertices convergent. In this case, there is no need 
for introducing a sharp momentum cutoff as in the local NJL model. Hence, the introduction of 
smoothing form factors in our expressions of an approximate local NJL model may be interpreted as 
a regularization procedure taking some effects of the originally non-local current-current interaction 
afterwards into account. 
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(1), which in addition to two independent coupling constants G%, G^ includes 
regularizing meson (diquark) form factors. Of course, it would be a very 
hard task to study the competition of %SB and CSC for arbitrary values of 
coupling constants G\, G2 and any form factors. Thus, in order to restrict this 
arbitrariness and to be able to compare our results (at least roughly) with 
other approaches, we find it convenient to investigate the phase structures 
of the model (1) at H = and H 7^ only for some fixed values of Gi, G2 
and some simple expressions for meson/diquark form factors (for simplicity, 
meson/diquark form factors are chosen to be equal). Let us choose the form 
factors f] 

A 4 A 4 
~ (A 2 +p 2 ) 2 ' 0n ~~ (A 2 +p 2 + #tfn) 2 ' (20) 

which have to be included in the energy spectra by a corresponding multi- 
plication of the a— , A— fields: 



E; = ^ + ^ : E r n = ^gHn+pl + <Pla\ |A'| - K\A\ (21) 
Note that with the choice of simple form factors ( |2D| ) our expression for the 



thermodynamic potential at H = formally coincides with the correspond- 
ing expression of Ref. |2Bj obtained for an NJL type model with instanton- 
induced four-fermion interactions. In particular, by a suitable choice of cou- 
pling constants Gi,G2, we will later "normalize" our phase portraits for 
H = to the curves of this paper in order to illustrate the influence of a 
nonvanishing chromomagnetic field. Q 

As a result, instead of fll^l) we shall deal with the following regularized 
potential V r (a, A): 



6 The application of the smooth meson form factors (20) leads in a natural way to a suppression 
of higher Landau levels, which is of particular use here. Hence, this regularization scheme is 
particularly suitable for the manifestation of the (chromo) magnetic catalysis effect of dynamical 
symmetry breaking. Indeed, the (chromo)magnetic catalysis effect and the underlying mechanism 
of dimensional reduction are closely related to the infrared dominance of the lowest Landau level 
with n = 0,|O|. 



7 It is necessary to underline that in our case the meson/diquark form factors (p0| ) mimick solutions 
of the BS-equation for some non-local four-fermion interaction arising from the one-gluon exchange 



approach to QCD. Contrary to this, the instanton-like form factor used in [28| has another physical 
nature. It appears as quark zero mode wave function in the presence of instantons |26||. 
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V r ((T,A) = V Q -2Nt(N c -2) J V 1 
where Vq = iV ( 



d 3 p gHN { °f ~ 



°^ oo 

/ Ean^dpa, (22) 



'« l 2d ^ G 2 



Ki=K + Tln 



' 1 + e -^ + M)V 1 + e -^-/i)' 



^2 = V(^ " J") 2 + 4|A|202 + ^ (£ r + ^2 + 4 | A |2 2 + 

+2Tln 



(23) 



^ + e -^-/i) 2 + 4|A^ / + e -^ + M )2 + 4 |A|2^ 

and £^ are given in ( |2T|) . Despite the A-modification, the expression ( p2|) 
contains yet UV-divergent integrals. However, as it was pointed out from the 
very beginning, we shall numerically study the subtracted effective potential, 
i.e. the quantity V r (a, A)-V^(0, 0), which has no divergences. 

In the next section the dependency of the global minimum point of the 
regularized potential ( |Z2| ) on the external parameters iif , /i, T will be investi- 
gated. 



IV. NUMERICAL DISCUSSIONS 



In the previous section we have chosen the form factors as in ( POf) in or- 
der to roughly normalize our numerical calculations at H ^ on the results 



obtained at H = in |28|] . Comparing the effective potential ( |Z2f) at gH = 



with the corresponding one from that paper (denoting their respective di- 
quark field and coupling constants by a tilde), we see that these quantities 
coincide if 2A = A, G\ = 2N C G\ and G2 = N c Gr2- Using further the nu- 
merical ratio of coupling constants from ref. [|28fl , we obtain in our case the 
following relation: 



G2 = t;G\. 



(24) 



Now, let us perform the numerical investigation of the global minimum point 
(GMP) of the potential ([22]) for form factors and values of coupling constants 
as given by ( |2D|) and (|2jj ), respectively. It was supposed earlier in some papers 



(see e.g. [|38jl ) that quantitative features of the color superconducting phase 



transition might indeed depend on the value of the form factor parameter 
A. So, in order to check the A-dependence of our results, we perform the 
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investigations for three different values of the cutoff, A =0.6 GeV, 0.8 GeV 
and 1 GeV. For each value of A, the corresponding value of G\ is selected from 
the requirement that the GMP of the function V# (cr, A)atT = fi = H = 
is at the point a = 0.4 GeV, A = in agreement with phenomenological 
results and [|28(| . Then, the value of G^ is fixed by the relation (|2|). This 
yields, for example, G X A 2 = 2N C 6A7 at A = 0.8 GeV, G X K 2 = 2iV c 6.16 at 
A = 1 GeV etc. 

Note, further that in the case of zero temperature we have 



Vi\ T =o = E r p + (fi-E r p )e{fi-E, 



V2\T=0 = - /i) 2 + 4|A|202 + ^( £ r + ^2 + 4 | A |2 2. 

In order to study the phase structure of the model using numerical meth- 
ods, the summation over n in fl22Q is limited by a maximum value n max = 



(2.5A) 2 /gH , where other terms of the series can be neglected due to their 
smallness. 

First of all, it should be remarked that, as in paper [^Bj at gH = 0, a 



mixed phase of the model was not found for H 7^ 0, i.e. for a wide range of 
parameters /i, H, T we did not find a global minimum point of the potential 
QZ2Q, at which a 7^ 0, A 7^ 0. The results of our numerical investigations of 
the GMP of V r (a, A) are graphically represented in the set of Figs. 1-6, where 
the notations I, II and III are used for the symmetric phase, for the phase 
with chiral symmetry breaking and for the CSC phase, respectively. For the 
points from region I the GMP of the potential lies at a = 0, A = 0. In region 
II we have a phase with broken chiral symmetry, corresponding to the GMP 
of the potential at a 7^ 0, A = 0. Finally, the color superconducting phase 
with the GMP of the potential at = 0, A 7^ 0, corresponds to the points 
from region III of these figures. 

In Fig. 1, one can see the phase portrait of the model in terms of /i, gH at 
T = for each of the above mentioned values of the cutoff A. The boundary 
between phases III and II is practically A-independent and represents a first 
order phase transition curve fi cr (gH). It is necessary to note also that for each 
value of A and a fixed value of gH there is a critical chemical potential jl c (H) : 
at which the GMP is transformed from a point of type III to a symmetric 
point of type I. However, this phase transition is remarkably A-dependent. 
Indeed, even in the simplest case with H = we have /2 C (0) = 1 GeV at A = 
0.6 GeV, /2 C (0) = 1.3 GeV at A = 0.8 GeV, /2 C (0) = 1.65 GeV at A = 1 GeV. 
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It is well-known from the one-gluon exchange approximation in QCD 
that CSC can exist even at enormously high values of the chemical potential 
/i ^ 10 8 MeV. So the above mentioned NJL framed transition from phase III 
to phase I looks, in the case T = 0, rather like an artefact of the regularization 
procedure. Q Thus, since this phase transition turns out to be unphysical, 
it is not shown in Fig. 1. The critical curves of this figure are obtained by 
interpolation in the most simple manner, i.e. by a second order polynomial, of 
numerical points lying at gH > 0.2 GeV 2 (for technical reasons, such points 
are explicitely shown in Fig. 2, rather than in Fig.l) Earlier, in the papers |I41 
the model (1) at G2 = and in the presence of an external magnetic field was 
considered. As shown there, for small values of the magnetic field strength the 
critical curves, as well as various thermodynamical and dynamical parameters 
of the system, demonstrate oscillating behaviour. In order to make more 
accurate interpolations and to become sure whether analogous oscillations 
appear in the present case (i.e., for the critical curves in Figs.l, 2) as well, 
or do not, one should make an enormous amount of numerical calculations, 
which proved to be rather difficult to accomplish. Due to this, we can make 
only a conjecture of an oscillating behavior of the critical curves judging from 
the positions of the points we have really calculated. Note further that in the 
region of low chromomagnetic fields, gH < 0.2 GeV 2 , we have extrapolated 
the critical curves to the known points at gH = 0. 

In Fig. 2, the phase portrait of the model in terms of fi,gH at T = 0.15 
GeV is presented for each value of parameter A =0.6 GeV, 0.8 GeV, and 1 
GeV. One can see that the boundary between II and III phases (a critical 
curve of a first order phase transition) only slowly changes with varying A. 
The second order phase transition from the CSC phase III to the symmetric 
phase I, which for finite T is now supposed to really exist, also has a weak 
A-dependence for gH < 0.2 GeV 2 . At greater values of gH the boundary 
between III and I phases has, however, a stronger A-dependence. The phase 
diagrams in the (T, fi) plane for gH = and gH = 0.4 GeV 2 are schematically 
represented in Fig. 3. The phase diagram in the (T, gH) plane for /i = 0.4 
GeV is represented in Fig. 4. For both figures we choose A = 0.8 GeV, for 



8 Note that in @ it 

was also claimed that in the NJL model at high enough the diquark 
condensate vanishes, which is the consequence of the regularization by a form factor. In this 
region, it might be necessary to use another approximation for the CSC investigation. 
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simplicity. It is necessary to point out that at gH = the numerical results 
of Figs. 1-4 coincide with those obtained in j28j at A = 0.8 GeV. Moreover, 
it should be emphasized that in all the above mentioned figures, a second 
order phase transition takes place at the boundary of the region I. At the 
boundary between regions II and III a first order phase transition takes place. 

Let us for a moment fix the value of the chemical potential and tempera- 
ture at varying values of gH. In this case, in the Figs. 1,2, one will have a 
straight line parallel to the gH axis. In particular, if fi = 0.4 GeV and T = 0, 
in Fig.l this line originates at gH = in the CSC phase III. At some (critical) 
value (gH) c pa 0.1 GeV 2 it crosses the line p cr (gH) and then, at yet greater 
values of gH, it passes through the phase II. Accordingly, at gH < (gH) c 
the GMP of the effective potential lies at the point (a = 0, A ^ 0), where 
A is equal to the diquark condensate in the true stable vacuum, whereas at 
gH > (gH) c the point (0, ceases to be a GMP. In this case it is only 

a local minimum point, so that A / corresponds to a metastable ground 
state of the system (for the stable ground state at gH > (gH) c the GMP is 
of the form a ^ 0, A = 0). Thus, the value (gH) c is the so-called evaporation 
point for the diquark condensate. 

In Fig. 5, the diquark condensate A is depicted as a function of gH 
for three values of the temperature at /i = 0.4 GeV and A = 0.8 GeV (in 
this Figure, due to problems with distinguishing closely positioned points for 
different cutoff values, we restricted ourselves to plotting curves for only one 
cutoff value A = 0.8 GeV.). Thick curves correspond to a stable diquark 
condensate (the point (0, A ^ 0) is the global minimum of the effective 
potential), and thin curves correspond to a quasistable diquark condensate 
(this is a local minimum). 

Here we should note that recent investigations yield the following value 
of the QCD gluon condensate at T = /i = 0: gH pa 0.6 GeV 2 J39j. It was 



shown in |40] that in the framework of a quark-meson model at ordinary nu- 
clear density po the gluon condensate decreases by no more than six percent, 
compared with its value at zero density. At densities 3po the value of (FF) 
decreases by fifteen percent. This means that for values of the chemical po- 
tential p < 1 GeV the gluon condensate is a slowly decreasing function vs 
/i. Taking in mind this fact, one can draw an important conclusion from our 
numerical analysis: At H = 0, T = and /i = 0.4 GeV there should exist the 
CSC phase (see [p3j|). However, if the condensate value gH pa 0.5 GeV 2 is 
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taken into account at fi = 0.4 GeV, then our model consideration concludes 
that the CSC phase does not exist at T = for such a large value of the 
gluon condensate (cf. Figs. 4, 5). f\ Assuming that our results would remain 
valid also for more realistic condensate fields, this would seemingly render it 
difficult to observe a CSC phase for T = at ji = 0.4 GeV. 

Notice, however, that our results change at finite temperature. First, if 
one takes T = 0.1 GeV and ji = 0.4 GeV, then at sufficiently small values 
of gH there is a symmetric phase of the theory, where both < qq > and 

< qq > condensates are zero (see Fig. 5). However, at the point gH « 0.1 
GeV 2 there is a second order phase transition to the CSC phase (here only 

< qq >^ 0), and at the point gH « 0.2 GeV 2 the external field destroys 
the CSC in favour of the chiral phase, where < qq >^ and < qq >= 0. A 
similar behaviour is observed at T = 0.15 GeV and ji = 0.4 GeV (see Fig. 5). 
Secondly, if one takes T = 0.15 GeV and \x > 0.6 GeV, then at sufficiently 
small values of gH there is a symmetric phase of the theory, where both 

< qq > and < qq > condensates are zero (see Fig. 2). However, at the 
point gH > 0.3 GeV 2 there is a second order phase transition to the CSC 
phase, where only < qq >^ 0. Hence, in some cases the gluon condensate 
can induce the CSC (the so-called chromomagnetic catalysis effect of CSC). 

This effect is further illustrated in Fig. 6, where the diquark condensate 
vs gH is depicted at ji = 0.8 GeV for two temperatures T = 0, T = 0.15 
GeV and two values of A = 0.8 GeV and A = 1 GeV. One can easily see that 
the value of A(gH) for T = 0.15 GeV is identically zero for gH < 0.44 GeV 2 
with A = 0.8 Gev, and for gH < 0.52 GeV 2 with A = 1 Gev, i.e., for each 
value of A, there exists a critical value of gH at which the nonzero diquark 
condensate is generated. 

Finally, let us make some additional remarks concerning the diquark con- 
densate at \i = 0.8 GeV. As it follows from our numerical analysis at T = 
(see Fig. 1, 6), for ji = 0.8 GeV, gH = 0, the GMP of the effective potential 
corresponds to the CSC phase with a stable diquark condensate A < 0.1 
GeV. However, assuming that the value of the gluon condensate gH « 0.4 
GeV 2 would hold for the above nonvanishing chemical potential, one would 
get a value of the diquark condensate A ^ 0.2 GeV, which is significantly 



'Recently, a similar prediction, namely that nonperturbative gluon fluctuations might be strong 



enough to destroy the CSC, was done in |4l|], but in a rather qualitative form. 
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larger in magnitude, than at gH = 0. It follows from Fig. 6 that the diquark 
condensate noticeably depends on A. One could suppose that this is due 
to the rather high value of the considered chemical potential, fi = 0.8 GeV. 
However, as our calculations show, a similar A-dependence of the diquark 
condensate occurs for smaller values of fi as well. These results confirm the 
conclusions of the paper [j3S[] obtained at H = that the value of the diquark 



condensate varies with A, if the form factor regularization is used in the NJL 
model. In contrast, the points of the boundary between II and III phases do 
not show a remarkable A-dependence (see Fig. 1,2). 

As a general conclusion, we see that taking into account an external chro- 
momagnetic field at least in the form as considered in the model above, 
might, in principle, lead to remarkable qualitative and quantitative changes 
in the picture of the diquark condensate formation, obtained in the frame- 
work of NJL models at H = 0. This concerns, in particular, the possibility 
of a transition to the CSC phase and diquark condensation at finite tem- 
perature. Clearly, a detailed quantitative discussion would, however, require 
to have additional information on the gluon condensate as a function of the 
temperature and chemical potential and on a possible /i-and giif-dependence 
of the cutoff parameter A. Q A further interesting generalization could be 
to extend this kind of approach to inhomogeneous background field config- 
urations [tefl. This concerns, in particular, the nonabelian condensate fields 



of the Stochastic Vacuum Model of QCD realizing Wilson's area law of con- 
finement |I3[ . 



V. SUMMARY AND CONCLUSIONS 

In the present paper the influence of different physical factors on the phase 
structure of the two-flavor NJL model (1) with two independent structures 
of four-quark interactions has been considered. This model is adequate for 
the description of the low energy physics of two-flavor QCD both in qq- and 
qq channels. In the papers [^7|- [|2T| it was shown that in QCD-motivated 



type of models (1) with H = the new color superconducting (CSC) phase 
can exist for moderate values of the chemical potential (baryon density). 



10 Such a dependence might, for example, arise, if one identifies the cutoff in the local 4-fermion 
model with the effective gluon mass in the gluon propagator. 
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As generalization of the "free field" NJL model of ref. [p9j| , we have, in 
particular, taken into account such nonperturbative feature of the real QCD 
vacuum as the nonzero gluon chromomagnetic condensate < F® v F aiiV >= 
2H 2 , which in the framework of a NJL model can be simulated by an external 
chromomagnetic field. The recent estimates give the following fixed value 
gH p hy S ~ 0.6 GeV 2 for the gluon condensate in QCD at T = 0. Despite 
of this fact, we considered it, however, useful to treat H as a free external 
parameter of the model. 

Since in the CSC phase the original SU C (3) symmetry of the theory is spon- 
taneously broken down to SU C (2), five color gauge bosons acquire masses. 
The corresponding external fields are expelled from the CSC phase (Meiss- 
ner effect). However, the other three "color isospin" gauge bosons remain 
massless, in accordance with the residual SU C (2) symmetry of the vacuum. 
Clearly, the corresponding external fields may then penetrate into the CSC 
phase. It is just the influence of these types of external chromomagnetic 
fields on the formation of CSC which was studied in our previous paper [j3T 



in the framework of a (2+l)-dimensional NJL model for vanishing chemical 
potential and temperature. There it was shown that chromomagnetic fields 
may induce the CSC phase transition. In the present paper, the chromo- 
magnetic generation of CSC has been studied in the framework of a (3+1)- 
dimensional NJL model for finite chemical potential and temperature. The 
vector-potential of the external chromomagnetic field was chosen to be of 
the Matinyan-Savvidy form (0) and lies in the algebra of the residual SU C (2) 
group, too. 

The coupling constants G\,G<i of our model (1) are considered as free 
independent parameters. In our numerical estimates, we found it convenient 
to use the relation (|2^ ) in order to "normalize" our calculations at H ^ 
on the known results at H = [|Z8j]. However, we hope that our qualitative 



conclusions remain also valid for values of G\ : G^ in some neighbourhood of 
The results of numerical investigations of the effective potential ( fZ2[) are 
presented in a set of Figs 1-6, where phase diagrams for the extended NJL 
model (1) in terms of /i, T, H as well as the behaviour of diquark condensates 
versus gH are shown. 

First of all we should note that the form factor regularization of the NJL 
model (1) is used throughout in the present paper. As shown in [|38|1 , in this 
case even at H = the features of CSC depend on the value of the form 
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factor parameter A. So, in order to clarify the corresponding situation at 
H 7^ 0, we have used three different values of A =0.6 GeV, 0.8 GeV, and 
1 GeV in the regularization scheme under consideration. One can see that 
phase transitions between the chirally broken and CSC phases have rather 
weak A-dependence (see Figs 1,2), whereas the diquark condensate values 
noticeably depend on A (see, e.g. Fig. 6). 

The main conclusion of our investigations is that the inclusion of an ex- 
ternal chromomagnetic field can significantly change the phase portrait, ob- 
tained at H = 0. Indeed, at H = T = the values of the chemical potential 
corresponding to the CSC phase approximately lie in the interval \x > 0.3 
GeV (see Fig. 1). If the external chromomagnetic field of the type (JTJ) is 
switched on at \i = 0.4 GeV, then at gH c « 0.1 GeV 2 there is a transition 
of the system from CSC to a phase, where only chiral symmetry is broken 
down. Thus, at gH phys « 0.6 GeV 2 , /i = 0.4 GeV and T = the CSC can 
not be observed at all. However, if T = and the chemical potential is fixed 
at \i = 0.8 GeV, then at least for all values < gH < 0.6 GeV 2 one can 
observe the CSC phase in which the diquark condensate A (if) is nonzero. 
It is worth mentioning that in this case the function A(H) is monotonically 
increasing (see Fig. 6) and the value of the diquark condensate at gH p h ys 
turns out to be significantly greater in magnitude than at vanishing H. Q 

Finally, one should note that at fi ^ 0, T ^ the external chromomag- 
netic field can induce the CSC phase transition. For example, at T = 0.15 
GeV and \i = 0.8 GeV there is a symmetric phase of the theory in which 
a = A = (both chiral and diquark condensates are zero) for sufficiently 
small values of gH (see Fig. 2). However, at some critical point gH c , a 
phase transition of the second order from the symmetric to the CSC phase 
is induced by the external chromomagnetic field. Notice that the CSC in- 
duction by some types of external chromomagnetic fields was observed in the 
framework of a NJL model at zero fi, T (see |Trj]- |33f). The present analysis 



shows that this effect takes place at some nonzero values of /i, T, too, which 
in principle could be important for heavy ion-ion collisions taking place at 



11 We roughly suppose throughout the present paper that at (fi, T) ^ the real gluon condensate 
is the same as at (/x, T) = 0. However, using a given fj,, T-dependency of the gluon condensate, it 
would be possible to extract physical informations about the CSC phase using our phase diagrams 
in Figs. 1-4 and plots of A(H) functions in Figs. 5, 6. 
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nonzero temperature. Clearly, a somewhat unpleasant feature of the above 
NJL approach, which requires some caution, is the A dependence of some of 
the results. Nevertheless, we believe that the above results are interesting 
and may serve as a starting point for further investigations of this issue. 
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APPENDIX: 



Let us sketch the calculation of the integral in eq. ([18]) of the text. Evi- 
dently, it can be rewritten contour integral 
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The integral over Cq is equal to the integral along the contour C (Fig. 8), 
and hence 

duo r da 2 (3 (3lu 
tg 



n 2 



1//3 2 



da 

(3a 

Y 



Co 2ni^ 2 a 2 + u 2 2 - 2 



2 / duj (3 



c 





1 


1 


) 1 


fc 


— ia 


u + ia 


/ 2ia 


(3a 


a=Q 


= 2 Inch 


n/3 


~2 


a=l//3 




2 ' 



0U> 



(A3) 



21 



REFERENCES 



[I] Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122, 345 (1961); 124, 246 
(1961). 

[2] D. Ebert and M.K. Volkov, Yad. Fiz. 36, 1265 (1982); Z. Phys. C 16, 
205 (1983); D. Ebert and H. Reinhardt, Nucl. Phys. B 271, 188 (1986); 

D. Ebert, H. Reinhardt and M.K. Volkov, Progr. Part. Nucl. Phys. 33, 1 
(1994). 

[3] D. Ebert, L. Kaschluhn and G. Kastelewicz, Phys. Lett. B 264, 420 
(1991). 

[4] U. Vogl, Z. Phys. A 337, 191 (1990); U. Vogl and W. Weise, Progr. Part. 

Nucl. Phys. 27, 195 (1991). 
[5] D. Ebert and L. Kaschluhn, Phys. Lett. B 297, 367 (1992); D. Ebert and 

T. Jurke, Phys. Rev. D 58, 034001 (1998). 
[6] H. Reinhardt, Phys. Lett. B 244, 316 (1990); A. Buck, R. Alkhofer and 

H. Reinhardt, Phys. Lett. B 286, 29 (1992). 

[7] N. Ishii, W. Benz and K. Yazaki, Phys. Lett. B 301, 165 (1993); C. 

Hanhart and S. Krewald, Phys. Lett. 344, 55 (1995). 
[8] J. Hiifner, S.P. Klevansky and P. Zhuang, Ann. Phys. 234, 225 (1994); 

I. N. Mishustin et al, Phys. Rev. C 62, 034901 (2000); |hep-ph/0010223| ; 
M. Hanauske, L.M. Satarov, I.N. Mishustin et al, Phys. Rev. D 64, 
043005 (2001). 

[9] S. Weinberg, Phys. Rev. D 13, 974 (1976); I.T. Dyatlov, JETP 82, 181 
(1996); see also the review by K. Yamawaki, in Proc. of "Dynamical 
Symmetry Breaking and Effective Field Theory", Cheju, Korea (1995) 
( |hep-ph/9603293| ). 

[10] T. Inagaki, T. Muta and S.D. Odintsov, Progr. Theor. Phys. Suppl. 127, 
93 (1997); A.S. Vshivtsev, MA. Vdovichenko and K.G. Klimenko, JETP 
87, 229 (1998); A. Goyal and M. Dahiya, J. Phys. G 27, 1827 (2001). 

[II] G.W. Semenoff LA. Shovkovy and L.C.R. Wijewardhana, Mod. Phys. 
Lett. A 13, 1143 (1998); W.V. Liu, Nucl. Phys. B 556, 563 (1999); E. 
Babaev, Phys. Lett. B 497, 323 (2001); Int. J. Mod. Phys. A 16, 1175 
(2001); V.Ch. Zhukovsky et al, Pisma Zh. Eksp. Teor. Fiz. 73, 137 (2001); 

E. J. Ferrer, V.P. Gusynin and V. de la Incera, jhep-ph/UlUlaua| . 

[12] S. Kawati and H. Miyata, Phys. Rev. D 23, 3010 (1981); V. Bernard, 
U.-G. Meissner and I. Zahed, Phys. Rev. D 36, 819 (1987); Chr.V. Chris- 
tov and K. Goeke, Acta Phys. Pol. B 22, 187 (1991); D. Ebert, Yu.L. 



22 



Kalinovsky, L. Miinchow and M.K. Volkov, Int. J. Mod. Phys. A 8, 1295 

(1993); A.S. Vshivtsev, V.Ch. Zhukovsky and K.G. Klimenko, JETP 

Lett. 64, 338 (1996); JETP 84, 1047 (1997). 
[13] S.P. Klevansky and R.H. Lemmer, Phys. Rev. D 39, 3478 (1989); H. 

Suganuma and T. Tatsumi, Ann. Phys. 208, 470 (1991); V.P. Gusynin, 

VA. Miransky and LA. Shovkovy, Phys. Lett. B 349, 477 (1995); T. 

Inagaki, S.D. Odintsov and Yu.I. Shil'nov, Int. J. Mod. Phys. A 14, 481 

(1999); E. Gorbar, Phys. Lett. B 491, 305 (2000). 
[14] D. Ebert et al, Phys. Rev. D 61, 025005 (2000); Phys. Atom. Nucl. 64, 

336 (2001); MA. Vdovichenko, A.S. Vshivtsev and K.G. Klimenko, Phys. 

Atom. Nucl. 63, 470 (2000); 
[15] K.G. Klimenko, Teor. Mat. Fiz. 89, 211 (1991); 90, 3 (1992); Z. Phys. 

C 54, 323 (1992); I.V. Krive and SA. Naftulin, Phys. Rev. D 46, 2737 

(1992). 

[16] V.P. Gusynin, VA. Miransky and LA. Shovkovy, Phys. Rev. Lett. 73, 
3499 (1994); Phys. Rev. D 52, 4747 (1996); Nucl. Phys. B 563, 361 

(1999) ; G.W. Semenoff, LA. Shovkovy and L.C.R. Wijewardhana, Phys. 
Rev. D 60, 105024 (1999). 

[17] VA. Miransky, Progr. Theor. Phys. Suppl. 123, 49 (1996); V.P. Gusynin, 
Ukrainian J. Phys. 45, 603 (2000). 

[18] D. Ebert and M.K. Volkov, Phys. Lett. B 272, 86 (1991); D. Ebert, 
Yu.L. Kalinovsky and M.K. Volkov, Phys. Lett. B 301, 231 (1993). 

[19] K.G. Klimenko, B.V. Magnitsky and A.S. Vshivtsev, Nuovo Cim. A 107, 
439 (1994); Theor. Math. Phys. 101, 1436 (1994); Phys. Atom. Nucl. 57, 
2171 (1994); LA. Shovkovy and V.M. Turkowski, Phys. Lett. B 367, 213 
(1996); A.S. Vshivtsev et al, Phys. Part. Nucl. 29, 523 (1998). 

[20] D. Ebert and V.Ch. Zhukovsky, Mod. Phys. Lett. A 12, 2567 (1997). 

[21] B.C. Barrois, Nucl. Phys. B 129, 390 (1977). 

[22] S.C. Frautschi, "Asymptotic freedom and color superconductivity in 
dense quark matter" , in Proceedings of the Workshop on Hadronic Matter 
at Extreme Energy Density, Ed., N. Cabibbo, Erice, Italy (1978). 

[23] D. Bailin and A. Love, Phys. Rept. 107, 325 (1984). 

[24] D.T. Son, Phys. Rev. D 59, 094019 (1999); D.K. Hong, Nucl. Phys. B 
582, 451 (2000); S.D.H. Hsu and M. Schwetz, Nucl. Phys. B 572, 211 

(2000) ; R.D. Pisarsky and D.H. Rischke, Phys. Rev. D 61, 074017 (2000); 
LA. Shovkovy and L.C.R. Wijewardhana, Phys. Lett. B 470, 189 (1999). 



23 



[25 
[26 



[27 

[28 
[29 

[30 

[31 
[32 

[33 
[34 

[35 



[36 

[37; 



[38 
[39 



[40 

[41 

[42 
[43 



K. Rajagopal and E. Shuster, Phys. Rev. D 62, 085007 (2000). 
R. Rapp, T. Schafer, E.A. Shuryak and M. Velkovsky, Phys. Rev. Lett. 
81, 53 (1998); Ann. Phys. 280, 35 (2000); M. Alford, K. Rajagopal and 
F. Wilczek, Phys. Lett. B 422, 247 (1998). 

M. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537, 443 (1999); 
K. Langfeld and M. Rho, Nucl. Phys. A 660, 475 (1999). 
J. Berges and K. Rajagopal. Nucl. Phys. B 538, 215 (1999). 
T.M. Schwarz, S.P. Klevansky and G. Papp, Phys. Rev. C 60, 055205 
(1999). 

M. Alford, |hep-ph/0102U47| ; B.O. Kerbikov, [hep-ph/0110T971 . 



D. Ebert, K.G. Klimenko and H. Toki, Phys. Rev. D 64, 014038 (2001). 
D. Ebert, K.G. Klimenko, H. Toki and V.Ch. Zhukovsky, |hep-th/01U5U49| 
(to be published in Progr. Theor. Phys.). 

D. Ebert, V.V. Khudyakov, K.G. Klimenko et al, |hep-ph/010818h . 



G.K. Savvidy, Phys. Lett. B 71, 133 (1977); S.G. Matinyan, G.K. 
Savvidy, Nucl. Phys. B 143, 539 (1978). 

N.K. Nielsen and P. Olesen, Nucl. Phys. B 144, 376 (1978); B 160 
380 (1979); N.K. Nielsen, M. Ninomya, Nucl. Phys. B 156, 1 (1979); J. 
Ambjorn, P. Olesen, Nucl. Phys. B 170, 60, 265 (1980). 
A. A. Sokolov, I.M. Ternov, V.Ch. Zhukovsky and A.V. Borisov, Gauge 
Fields [in Russian], Moscow University Press, Moscow, 1986. 
R.T. Cahill and CD. Roberts, Phys. Rev. D 32, 2419 (1985); R.T. 
Cahill, CD. Roberts and J. Praschifka, Phys. Rev. D 36, 2804 (1987); J. 
Praschifka, R.T. Cahill and CD. Roberts, Int. J. Mod. Phys. A 4, 4929 
(1989); D. Ebert et al, Nucl. Phys. A 552, 523 (1993). 
M. Sadzikowski, Mod. Phys. Lett. A 16, 1129 (2001). 
S. Narison, Phys. Lett. B 387, 162 (1996); M. D'Elia, A. Di Giacomo, E. 
Meggiolaro, Phys. Lett. B 408, 315 (1997); A.I. Alekseev and B.A. Arbu- 
zov, Mod. Phys. Lett. A 13, 1747 (1998); A.I. Alekseev, |hep-ph/98082TJ5 



K. Saito, K. Tsushima and A.W. Thomas, Mod. Phys. Lett. A 13, 769 

(1998) . 

N.O. Agasian, B.O. Kerbikov and V.I. Shevchenko, Phys. Rept. 320, 131 

(1999) . 

D.G. Caldi and S. Vafaeisefat, Phys. Lett. B 287, 185 (1992). 
H.G. Dosch Phys. Lett. B 190, 177 (1987); H.G. Dosch and Yu. A. 
Simonov, Phys. Lett. B 205, 339 (1988); Yu.A. Simonov, Nucl. Phys. B 



24 



324, 67 (1989); H.G. Dosch, Prog. Part. Nucl. Phys. 33, 121 (1994). 



25 



Figure captions 



Fig. 1. Phase portrait of the model in terms of variables (fi,gH) at 
T = for three values of the cutoff parameter A = 0.6, 0.8, 1 GeV. Regions 
II, III describe the phase with broken chiral symmetry (a / 0,A = 0) and 
the color superconducting phase (a = 0, A ^ 0), respectively. 

Fig. 2. Phase portrait of the model in terms of variables (/j,,gH) at 
T = 0.15 GeV for three cutoff values A = 0.6,0.8,1 GeV. The included 
phase I is the symmetrical one {a = 0, A = 0). 

Fig. 3. Phase portrait of the model in terms of variables (T, li) at gH = 
and gH = 0.4 GeV 2 for A = 0.8 GeV. 

Fig. 4. Phase portrait of the model in terms of variables (gH, T) at 
fj, = 0.4 GeV and A = 0.8 GeV . 

Fig. 5. Diquark condensate as a function of an external chromomagnetic 
field for ji = 0.4 GeV, A = 0.8 GeV and three values of the temperature, 
T = 0, T = 0.1 and T = 0.15 GeV. 

Fig. 6. Diquark condensate as a function of an external chromomagnetic 
field for ji = 0.8 GeV and two values of the temperature, T = 0, and T = 0.15 
GeV and for two values of the cutoff parameter A = 0.8, 1 GeV. 

Fig. 7. Integration path Cq used in eq. (Al) of Appendix. 

Fig. 8. Integration path C used in eq. (A3) of Appendix. 
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